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ABSTRACT 

We introduce a simple generalization of the basic holographic superconductor model in 
which the spontaneous breaking of a global U{1) symmetry occurs via the Stiickelberg 
mechanism. This more general setting allows tuning features such as the order of the 
transition. The physical vacuum of the condensed phase and the order of the transition 
are determined by a detailed analysis of the free energy of the system. For first order 
transitions, we identify a metastable phase above the critical temperature. In this case, 
the conductivity shows additional poles, thus suggesting that the condensate has internal 
structure. We comment on the possibility of obtaining second order phase transitions with 
non mean-field critical exponents. 
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1 Introduction 

The AdS/CFT correspondence [U [2], [3] is a powerful tool for dealing with strongly cou- 
pled CFTs through a weakly coupled dual gravitational description^ Conformal field 
theories are relevant in condensed matter physics near quantum critical points. Hence, 
the AdS/CFT has the potential to provide a handle on the dynamics of quantum critical 
points which cannot be attained by other analytical methods (see the excellent reviews 
O [HI [7] and references therein for a detailed introduction). More generally, close to second 
order phase transitions, correlation lengths diverge and systems are well approximated by 
a continuum scale invariant theory. In this paper, we take a more general approach and 
consider CFTs even when the previous conditions are not met, for example for first or- 
der phase transitions. In this case, the CFTs should be regarded as examples of strongly 
coupled gauge theories with fairly simple gravity duals. 

More specifically, we investigate strongly coupled CFT^_i's that exhibit spontaneous 
symmetry breaking below some critical temperature along the lines of |H1 [H]- Most of 
the details of this CFT^.i, other than the existence of a global U{1) symmetry which is 
spontaneously broken by the Stiickelberg mechanism are largely irrelevant. We will 

^Non-conformal extensions of the correspondence are also known, with [4J the prototypical example. 
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just assume that the CFT(i_i dynamics can give rise to spontaneous symmetry breaking of 
the U{1) symmetry controlled by an order parameter given by some scalar operator O of 
conformal dimension in the CFT^.i. Then, as a function of temperature, the system 
will be either in an ordered (broken symmetry) phase at low temperatures or in a disordered 
(symmetric) phase at high temperatures. We will assume the CFTrf_i has gravity dual, in 
terms of which we compute the properties of the system. The study of this phase transition, 
along with the transport properties of the system, is usually beyond standard analytic field 
theory methods but still accessible by using the AdS/CFT correspondence. 

The central observation of this paper is that, under the very generic premises demanded 
on the field theory side, namely a CFTd_i with a global U{1) symmetry which can be 
spontaneously broken, the dual gravitational description can be considerably more general 
than the set-up considered in [HI [9]. Taking an effective field theory approach, we will 
propose a class of gravity duals with the right ingredients to model the considered CFTs. 
We expect that all these gravitational duals do indeed yield a phase transition in a similar 
manner to [HI |9]0 However, the more general setting provided by the Stiickelberg model for 
spontaneous symmetry breaking permits tuning certain features of the phase transition, 
such as its order. Thus our results provide an effective description of different phase 
transitions in strongly interacting systems induced by the spontaneous breaking of a global 
f/(l) symmetry]^ 

The outline of this paper is as follows. In section 2 we introduce a gravity dual of 
the Stiickelberg model for spontaneous symmetry breaking which depends on a general 
function JF of the scalar field. We then find the equation of motion of the gravity dual for 
the simplest choices of JF, namely, monomials of degree n of the field. Finally we provide 
qualitative arguments, in very much the same spirit as in [S], that a phase transition should 
occur. In section 3 we prove numerically that the transition indeed exists and that its order 
depends on n. For n = 2 the transition is second order and for n > 2 it is first order. In 
both cases we find that there are actually many solutions to the equations of motion. The 
most stable -and therefore physical- solutions are determined in section 4 by analyzing the 
free energy. In the case of first order phase transitions we have identified the expected 
metastable regime around the critical temperature. The transport properties of some of 
these CFTs are studied in section 5. For first order phase transitions (n > 2) we find extra 
poles in the conductivity. This hints the condensate in the ordered phase has internal 
structure. In section 6, we investigate more general choices of !F and discuss the possibility 
of non mean-field critical exponents. We end with some conclusions in section 7. 

^Other interesting realizations of holographic superconductors have been introduced in [iTJ [12l [13l [Mj 

nsjiiniiniiTH]. 

•^We are considering that the U{1) symmetry will eventually be gauged. Strictly speaking, a theory with 
a spontaneously broken global U{1) symmetry corresponds to a superfluid. The holographic description 
of such system was studied in JT^ I^Uj • 
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2 Stiickelberg model for spontaneous symmetry break- 
ing 

We are interested in strongly coupled CFT(i_i's which enjoy a global U{1) symmetry. We 
suppose that the strong coupling dynamics of the CFT can lead to spontaneous symmetry 
breaking of such a symmetry. This symmetry breaking is controlled by some order 

parameter. We will assume it to be some scalar operator O of conformal dimension in 
the CFTrf_i. Since we assume this CFTd_i has a gravity dual, on very general grounds, the 
global U (1) symmetry gets mapped to a bulk U{1) gauge symmetry. Likewise, the operator 
O translates to a bulk scalar field \1/ with the appropriate mass given by A^. Therefore, the 
minimal requirements to holographically model the systems of interest is a f/(l) gauge field 
and a scalar; both coupled to gravity and living in an (asymptotically) AdS^ background. 
Since we are interested in the finite temperature case, the background will actually be a 
black hole in AdS^. Spontaneous breaking of the boundary U{1) symmetry at some critical 
temperature translates in the holographic dual into condensation of the scalar. This leads 
to spontaneously breaking of the bulk gauge symmetry through the Higgs mechanism, 
which is dual to the spontaneous breaking of the boundary global symmetry [2T1 [22] . 

The condensation of the scalar actually changes the background from the AdS-black 
hole in the uncondensed phase to a hairy black hole in the condensed phase. The following 
minimally coupled scalar [H [9] is one of the simplest models with this property. 



S 



f ( F'^ \D \I/P ^ 

J v^{R-^--^-^-^-^^-vm)}, (1) 

where = — iA^. We can innocuously re-write this model in a Stiickelberg form by 
re-writing the charged scalar field \1/ as \&e*^: 

r ( P2 2^2 ^2 

S = j ^9{R-^--^'-^-^-\{dp-Af-V{^>)] (2) 

with \1/ and p real. The gauge symmetry becomes A — > A + da together with p p + a. So 
far all we have done is rewriting the model. Nevertheless, it is straightforward to generalize 
the model preserving gauge invariance. The generalized action reads 

r ( P2 n^2 ^2 ^2 ^ . 

s = j v-9{--i--^-^-\m)\{dp~Af-vm] (3) 

where JF is a function of ^ . We take its absolute value to ensure positivity of the kinetic 
term for We will refer to this model as the Stiickelberg holographic superconductor. Our 
generalization of the basic holographic superconductor takes a rather compact form. 



^One could in principle be worried about the non-analyticity of such an action due to the absolute 
value. In all the examples we study in the paper, the solutions remain well inside the region oi J- > 
for all radial positions. We can always think about modifying (if necessary) such that the solutions are 
preserved but \J- \ remains non-zero and analytic for all 
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When re-writing ([l) in the form (|2|, it is important to remember that ^ is the absolute 
value of a complex field and hence positive definite. Keeping the definition of our model as 
broad as possible, ^' in (j3| takes in general any real value. Nevertheless, as it happens for 
(|2|, it is sometimes natural to restrict ^ to be positive if the model admits a reformulation 
in terms of a complex field. We will exploit this fact later. The results in [HI [9] strongly 
suggest the detailed form of the potential V are largely irrelevant as far as the condensation 
is concerned. Thus, for the sake of simplicity, we will assume ^ = from now on. 



2.1 Choice of 

In order to determine the action, we need to specify the function JF in (|3|. If JF is an 
analytic function of \E', it admits a Taylor expansion. In what follows, we will focus mostly 
on some simple models in which JF is just a monomial of some degree n 

T ~ (4) 

This choice is a good starting point to initiate the analysis of Stiickelberg superconductors. 
In section 6, we will explore more general choices of T . 

To be more specific, the theories we consider contain gravity with a negative cosmolog- 
ical constant A plus a real scalar field a real pseudoscalar field p and a f/(l) gauge field 
whose coupling constant we call e. The action reads, 

Sn = j rf"+^v^{i?-A+ (5) 

The scale M, required by dimensional analysis, will be set to the unityj^ 

This theory has the aforementioned local gauge invariance under which the p and Am 
fields transform as 

p-^p + a{x^) , Am Am + dMa{x^) . (6) 

This symmetry is broken if \^ develops a VEV. In flat space this would never occur without 
a Higgs-like potential for However in a non-trivial background, as in [8],|9], the presence 
of a chemical potential or a charge density might induce the condensation of \E'. Indeed 
upon a trivial numerical re-scaling, the n = 2 theory is just the one considered in [HI |9], 
where the existence of a condensed phase was explicitly shown. 

^Eventually these phenomenological models might come from a UV completion such as string theory. 
It is then natural to assume that M ~ L~^, being L the AdS radius. Since we will work in units in which 
L = 1, this provides justification to set Af = 1. 
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2.2 Equations of motion 

We are now ready to derive the equations of motion for the different fields in (|5]). We start 
with the gravitational field. The Einstein equations are sourced by the matter stress-energy 
tensor. The latter is given by 



-MN 



Qmn 
4 



8e2 



+ Qmn 



m ~ . 



+ 9MN ^ + 9MN ^" (^AP - Aa^ 



2e2 



After a re-scaling \^ = e\E', 

Tmn 1 / FabF"^^ 



■MN 



9mn 



+ 9mn 



Oa-^O^-^ 



m 



N 



(7) 



9mn ^ * + 9mn i^AP - Aa) 



FmbF^ n Om-^On^ 



^^{dMP-AM){dNP-A 



(8) 



In this form it is clear that in the probe limit (e oo) gravity is decoupled from the rest 
of fields. For the sake of simplicity we will work in this limit in the rest of the paper. 

Since there is a negative cosmological constant, we will consider AdS^+i black holes, 
whose geometry is given by 



ds' = -f{r)de + j^^+r'dxl_, , /(r) 



M 



d-2 



r'd—2 



1^2 J.U 

The dynamics of the other fields in this background is described by the action. 



(9) 



Sn = j C^'+'v^{ 



FmnF^^^ dM^d^^^ m2 
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(10) 



The equations of motion of these fields are given by the Lagrange equations coming 

(11) 



from (10) in the above AdS^+i background: 



and 



dsiyg"^'' (s^p-A^)) = 0, 







(12) 
(13) 



We just keep nllv and (13) since (12) can be obtained by acting with ds on (11). 



We use the gauge freedom to fix p = 0. In addition, we consider that only the time 
component Aq of the gauge field, which we will call $, is turned on. Since we are interested 
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in hair-like solutions we will also assume that both and $ are only functions of r. Then, 
particularizing for the background of interest, the relevant equations become, 

d — 1 

$" + $' _ 2 — $ = ; (14) 

r J 



^ + ( + ^)^'-— ^ + n— — $2 = ; (15) 



and 

/ P 

where the prime denotes derivative with respect to r. 

For n > 1 it is clear that a solution of the above equations with the required boundary 

conditions is \E' = and $ = /x This is just the normal (symmetric) phase. We will 

later show that for sufficiently low temperatures (~ r^j), a condensed solution representing 
the symmetry breaking phase exists. 

It is convenient to define the dimensionless coordinate z = where th is the horizon 
radius 

4 = L^M^-^ . (16) 

In these coordinates the horizon sits at z = 1 while the boundary of the AdS space is at 
z = 0. The equations of motion are given by, 

^ z z^l-z^) ' ^ ' 

and 

where now prime denotes derivative with respect to z. 



2.3 Boundary conditions 



In order to solve the equations above, we have to fix the behavior of the field at the horizon 
and the boundary. Close to the boundary the scalar field has free field behavior. Therefore, 



I^Z 



v]/_r 



H 



X- 



d±^d'^ + Am? 



L2 



(19) 



Since we are interested in spontaneous symmetry breaking, we focus on normalizable 
modes, such that their coefficient is proportional to the VEV of the dual operator 0± 



at the boundary. In order to interpret the solution to (17) and (18) as hair for the black 
hole the fields must remain regular everywhere. Normalizability of $ dt requires that $ 
vanishes at the horizon. Furthermore, regularity of the solution requires ^ to be finite at 
the horizon 

^1 = a + h{l- z) + --- (20) 



and to decay as (19) at the boundary. 
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Likewise, close to the boundary, $ behaves as 

<^-H-^,z'-\ (21) 

where is a chemical potential and p a charge density. We will later solve the equa- 
tions numerically and clarify whether our model has an ordered phase for sufficiently low 
temperature. In the next section we provide heuristic arguments that this is the case. 



2.4 Qualitative arguments for the existence of the transition 

We note that the equation of motion for ^ is the same as the one we would obtain from 
an effective Lagrangian 

Seff = j V9 (9m^5^^ + m2vl/2 + Ae// vl>-) , (22) 

where the effective coupling is 



(23) 



H 



Thus, we have effectively a theory for a scalar field with a given potential of order n (let 
us assume n > 1). A crucial observation is that this coupling is negative. These types of 
unbounded potentials are known to give rise to hairy black holes (see e.g. [231. 12 ^ 125]). This 
is our first indication that a condensed phase might actually exist for any n. Nevertheless, 
this is not sufficient to argue that a phase transition indeed exists. In this situation, the 
hair is in some sense more a property of the boundary than of the horizon. However, this 



is not the case for the potential above, since the coupling (23) depends on properties of 
the horizon through $2 (which depends on the chemical potential and the density), th and 
the extra redshift factors. Then, it is not surprising that, as we will later see, the potential 
is sufficiently strong to trigger the existence of hair, at least within a certain window of 
parameters. 

Let us now consider the \1/ equation of motion, which can be rewritten as, 

^" + ^ - J- im^ + ^ ^"-2 \ ^ 0. 24 

z{l-z'^) z^{l-z'^)\ rj^{l-z'^) J ^ ' 

Close to the boundary, the potential-like term is suppressed with respect to the mass term. 
This explains why the existence of a condensed phase is not a priori guaranteed, since 
the potential effectively does not extend all the way to the boundary. On the contrary, 
any non-vanishing expectation value for the scalar field at the horizon would lead to an 
instability in the electric field, in the same spirit of [U |9]. 

The value of n does indeed affect the details of the condensation. The reason is that 
\1/ goes from \l/o in the horizon to at the boundary. Since \1/ appears raised to the 
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power n — 2 in (24), it is natural to expect that the larger n is, the more sudden this 
contribution kicks in. From this heuristic argument, we expect that the phase transition 
becomes sharper as n increases (in fact, we will later see that it becomes first order for 
n > 2). Another consequence of this argument is that the critical temperatures should 
decrease for larger n. 



3 Numerical evidence for a phase transition 

In this section we use the shooting method to numerically solve the equations of motion 



(26) and (27) with the boundary conditions discussed in the previous section. In order 
to compare with p], we work in the canonical ensemble. In this ensemble, the density p 
is kept fixed and the chemical potential is determined by the solution of the equations of 
motion. Moreover, for concreteness, we assume d = 3 and mL^ = —2. We do not expect 
that other choices will qualitatively modify our results. For the above values of d and m 
the behavior of the scalar field at the boundary is 

^ = — 2 + ^2^ 25 

We note that the mass, though tachyonic, is above the Breitenlohner-Freedman (BF) bound 
in AdS4|^ In addition, both modes and ^2 are in this case normalizable. We define 
{Oi) = \E'j. In order to have a well defined solution, we have to choose a quantization in 
which only one of the {Oi) is non- vanishing. The analysis of the free energy in the next 
section will confirm this point. 

For the numerical calculation we absorb r^/ in $ as <y9 = ^ and define X = f • This 
choice is numerically more efficient since the equations do not depend explicitly on the 
temperature and the \1/ expansion at the boundary starts with a constant instead of with 
a z suppression. In terms of these new fields, the equations of motion read 

97-2 n-2 

V'"-^r3^X> = 0, (26) 

and 

^ ~ Z{l-Z^) ^ ~ (1-^3) ^ ^ (1-^3)2 X = . (7) 

In what follows, we focus on the cases n = 2 and 3. We have also investigated n > 3 
and observed that the main qualitative features are not changed. 

The condensed phase, where (Oi) 7^ 0, only exists below some temperature, which we 
denote Tq. The critical temperature T^, is the one at which the free energy of the condensed 
phase becomes smaller than the one for the un condensed phase. Generically, Tc < Tq. For 
second order phase transitions, Tc = Tq. These issues are discussed in detail in section |4| 



^In fact, as discussed in |27j . the mass term plays no important role in the existence of condensation for 
the scalar. The value of the mass affects various characteristics of the resulting superconductor, though. 
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In figure [T] we present results for tlie Oi condensate normalized by the critical tem- 
perature for n = 2,3. Figure |2] shows a similar plot for the normalized O2 condensate. 




Tc 



Figure 1: Normalized condensate {Oi) as a function of temperature at fixed p = — 1. 
The blue and magenta curves correspond to n = 2 and n = 3, respectively. For n = 2, 
Tc ~ 0.268 while for the n = 3, ~ 0.187. 

Similarly to what happens for = 2, it is possible to derive the n = 3 model from one 
with a complex scalar such that \1/ is constrained to be positivejj 




2 ^ 



0.5 0.6 0.7 0.8 0.9 1.0 

T 
Tc 

Figure 2: Normalized condensate ^y {O2) as a function of temperature at fixed p = —1. 
The blue and magenta curves correspond to n = 2 and n = 3, respectively. For n = 2, 
Tc ~ 0.141; while for n = 3, Te ~ 0.113. 

Figures [1] and [2] indicate that the transition is first order for n = 3 (preliminary analysis 
indicate that the same behavior occurs for n > 3). In other words, the condensate remains 

^If we allow to take any real value, both the n = 2 and 3 models have a 5* ^ — ^' symmetry that 
allows us to focus only on positive values. 
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finite at the temperature Tc at which the free energies of the condensed and uncondensed 
phases become equal. In section 4, we investigate this issue in greater detail. Indeed, in 
figures [T] and [2] the n = 3 curve is the upper half of the full solution. In figure [3} we show 
the full solution for the Oi condensate (a very similar plot holds for O2). Figure [s] shows 
that there are two possible values of the condensate for each temperature]^ Since the upper 
half of the curve has the largest condensate, it is natural to assume (we confirm this in 
section 4 by studying the free energy) that the curves in figures [l] and [2] correspond to the 
physical solutions. 




o' ^ 1^-— - I < -- 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

T 



Tc 

Figure 3: Normalized condensate (Oi) as a function of temperature at fixed p = —1 for 
n = 3. The thicker curve indicates the physical piece of the full solution. 

Finally we note that according to figures [l] and [2] (and similar ones for n > 3, that we 
do not exhibit here) the value of the condensate increases with n. This suggests that the 
attractive interaction that leads to the formation of the condensate becomes stronger as n 
is increased. 



3.1 The n = 2 model: a phase space for the system 



In this section, we investigate in more detail the solutions to (|26j) and (27) for n = 2 
and either Oi or O2 equals to zero. We will see that, for sufficiently low temperatures, 
these equations have multiple solutions or "branches" n An interesting way of visualizing 



multiple branches is to momentarily allow Oi 7^ and search for solutions to (26) and 



(27), still at fixed p. As we will see in the next section, one of the condensates must vanish 



in order to have a critical point of the free energy, which is a necessary condition for the 



^ A similar behavior was first observed in pjjj , in the case of holographic superfluids in the presence of 
a non-zero fluid velocity. 

^The analysis of the free energy of these solutions (see section |4| indicates the system condenses into 
the branch that shows up at the highest temperature, which also corresponds to the one with the largest 
condensate. 



11 



stability of the solution. Thus, out of the resulting curve in the {Oi,02) plane, we can 
focus just on the intersections with the axes, at which one of the Oi vanishes. We will later 
determine which of these branches is indeed stable, i.e. a minimum of the free energy. In 
some sense, the (Oi,02) plane can be regarded as a sort of phase space for the system, 
which contains interesting information. 

As we said, we now focus on = 2. Given that for Oi is bigger than for O21 we 
can identify three different regions depending on the temperature: 

Region I: Tf ^ < T Region II: Tf ^ < T < Tf ^ Region III: T < Tf ^ 
Region I 

In this region, T is bigger than any of the two critical temperatures. Therefore, we should 
expect no cuts with the axes. The resulting plot, which satisfies these expectations, can be 
seen on the left of figure |4| This behavior indicates that the curves in figures 1 and 2 are 
the first ones to appear as we decrease the temperature and hence determine the critical 
temperature T^. For T > T^. the system is in the non-condensed phase. 

Region II 

In this region, T is smaller than the critical temperature for Oi, but larger than for 02- 
We therefore expect a cut with the horizontal axis but no cut with the vertical 
shown on the right of figure |4j 




012345 012345 
<0i> <0i> 

r = 0.30 T=0.20 
Figure 4: (Oi, O2) plane at T = 0.30 (Region I) and at T = 0.20 (Region II). 
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Region III 

In this region, T is smaller than both T^'^ and T^^ . We then expect cuts with both axes. 
In fact, for sufficiently low temperatures, the spiral circles the origin multiple times. As 
the temperature is decreased, new turns emanate from the origin. As a result, for very 
low temperatures, there are multiple intersections with the axes (i.e. branches). Figure [s] 
shows spirals for different temperatures. 



<02> 




-0.2 0.2 0.4 0.6 

<0i> 

r = 0.050 

Figure 5: (Oi, O2) plane for T = 0.075, 0.065 and 0.050, well inside Region III. 

Figure |6] shows the corresponding O2 = ^ branches as a function of temperature. In 
section 4, we will show that the new branches have larger free energy than the ones con- 
sidered in figures 1 and 2, corresponding to the outermost cuts. Furthermore, while all 
the branches are critical points of the free energy (i.e. its derivative with respect to the 
vanishing Oi is zero), only the outermost branches correspond to local minima. 

As we have explained, the outermost solution corresponds to the one with the largest 
Tc. This suggests that it should be the one preferred by the system. This picture will be 
confirmed by the study of the system free energy. 

It would be interesting to perform the same analysis for general choices of the function 

and explore whether a behavior similar to the one of the n = 2 model holds in other 
examples. 
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Figure 6: |(9i| as a function of the temperature for the different branches. Notice that, 
in order to simphfy the comparison with figure |5] we have not normahzed the condensate 
by the temperature as in previous plots. The dashed hues indicate the three temperatures 
considered in figure [5] 

4 The free energy 

In this section we analyze the free energy for n = 2,3. We first review how to compute the 
free energy in models with gravity duals. 



4.1 General considerations 

The free energy of a system is given by its on-shell action. In the models we are interested 
in, this simple procedure involves a few subtleties that we address following [19] . 

Let us start with a little digression about theories in spaces with boundary. In this case, 
there is a boundary term in addition to the equations of motion in the bulk. So far we have 
not been careful about this contribution, assuming that the necessary boundary counter- 
terms were suitably added. To have a well-defined variational problem, not only the bulk 
equations of motion must be satisfied, but also the boundary term must vanish. Usually, 
for a theory in a space without boundary, this is achieved by the physically reasonable 
demand that fields vanish at infinity. However if the space where the theory lives has a 
boundary we have to tackle the interplay of the boundary conditions and the boundary 
term. In our case, the boundary term is 



/ \y^(FrtS^-g'-'"dr^S^)] (28 

J Boundary ^ ^ / ^ Boundary 



' Boundary 

We can take Boundary = 0. Since 5$ = 5/i at the boundary, this choice corresponds 
to the grand canonical ensemble, in which /x is fixed. If, in addition, we choose boundary 
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conditions such that = 0, the second term also vanishes. The boundary value of \l/ is 
fixed and, consequently, it corresponds to the choice Oi quantization. 

Under these assumptions we have a well-defined variational problem. Nevertheless, the 
action is still divergent. In order to regularize it, we introduce a cut-off before the 
boundary (r oo), which we will eventually send to infinity. One can verify that 



^onsheii - 2 + 2L + 2L 2 I zHl-z^) ^ ' 

where the superscript (1) stands for the choice of Oi quantization. 

In order to carry out a holographic renormalization of the free energy, we must add the 
boundary counter-term, 

{ASBoundaryY'^ = -I [ [Vl Boundary ^ 

J Boundary 

where 7 is the induced metric at the boundary. 

The renormalized free energy W'^^^ = -S^^l ^f^^^ - {ASBoundaryY^^ is given by 



2 2 2 Jo 1-z^ ^ ' 

Once again, the superscript indicates that we have chosen the Oi quantization. 

As explained in j2T], Oi and O2 are canonically conjugated variables. We can see this by 
fixing 5\l/2 rather than in our general variational problem. In order to do this, we note 



that, in the grand canonical ensemble, the boundary action (28 ) reads [—r% S'i/] Boundary 
Thus, it can be cancelled by the variation of 



S 



Boundary 



Boundary 



(32) 

Boundary 



Therefore, by considering SBuik + SBoundary, we have a well-defined variational problem 
provided \l/2 is fixed at the boundary. It is easy to check that the on-shell action is 
again divergent. As before, it can be renormalized by adding the appropriate counterterm 
(A5)(2) = -(AS)(i), resulting in 



,VC. + ^ + li^^ll /' f^W. (33) 

2 2 2 Jo 1-^3 

Then, we can think about ly^^) ^j-^g Legendre transformation of W^^\ In particular. 

Thus, the local extrema of W^'^^ sit at vanishing O2 and, conversely, the extrema of VT^^) g^Ye 
at Oi = 0. This justifies our previous assumption that only one of the is non-vanishing 
for physical solutions. 
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Similarly we can keep p, instead of /i, fixed at the boundary, i.e. work in the canonical 
ensemble [19]. In what follows, we focus on the grand canonical ensemble with /i = — 1. 
Moreover, we mostly present results for the Oi quantization. 



4.2 Free energy for the n = 2 model 

We are now ready to compute the free energy as a function of temperature by evaluating 



(31) with fields given by the solution to the equations of motion. An interesting quantity 



is the difference in free energies between the condensed and uncondensed phases. 



Condensed Uncondensed' 



(35) 



where Wjjl^^^^^^^^ is the free energy for = \E'2 = 0. To get a flavor of the behavior of 
the free energy, we can extend our previous plots by adding it as third dimension to the 
spirals in the {Oi,02) plane, as shown in flgure [7j The free energy typically increases as 
we approach the origin, i.e. the uncondensed phase. A more thorough numerical analysis 
shows that the external branch (in the case of AW^^\ the branch associated with the 
outermost cut between the spiral and the 02 = axis) is a local minimum of the free 
energy. On the other hand, internal branches are critical points but not extrema. In this 
way, we conflrm the previous statement that only the external branch corresponds to a 
condensed phase. The local minimum becomes more manifest as we approach the critical 
temperature. 

<02> 




-0.15 -0.1 -0.05 

<0i> 



0.05 




Figure 7: Normalized AW*^^) as a function of {Oi, O2) for n = 2, fi = —1 and T = 0.015. 

Figure |8] shows the free energy as a function of the condensate around the critical 
temperature. We observe the characteristic behavior of a second order phase transition. 
Above Tc, the only minimum is at the origin. Below Tc, the origin becomes a maximum and 
a minimum develops at non-zero condensate, which smoothly approaches zero as T — > T^. 
The minimum in the Oi direction occurs at {O2) = 0. 

Finally, flgure [9] shows the smooth convergence of the free energies of the condensed 
and uncondensed phases as the temperature approaches Tc 
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Figure 8: Normalized AVT^^^ as a function of {Oi) for n 
0.20,0.25,0.30,0.35. 
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Figure 9: Normalized difference in free energy between the uncondensed and condensed 
phases as a function of the temperature. 



4.3 Free energy for the n = 3 model 

We can use our results on the free energy to get a better understanding of the phase 
transition in the ra = 3 case. An important difference between tt, = 3 and = 2 is the 
existence of a lower branch in figure |3} Figure 10 shows the free energy for both branches. 
The branch with the largest condensate has less free energy. Furthermore, it is a local 
minimum of the free energy, while the other branch is a local maximum. We conclude 
that, as we previously assumed, the system prefers to condense in the branch with the 
largest condensate. 



Figure 11 shows the free energy as a function of the condensate for various tempera- 
tures. Below some temperature Tq, a metastable minimum develops at a finite value of the 
condensate. As we lower the temperature, we reach the critical temperature Tc, at which 
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0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 

T 

Figure 10: Normalized AVT*^^^ for the two branches as a function of the temperature. 

this minimum becomes degenerate with the uncondensed one. Since, at Tc, this minimum 
is at a finite distance from the origin, the transition is first order. For T < Tc, the system 
is in the condensed phase. 




0.5 1 1.5 2 2.5 

<0i> 



Figure 11: Normahzed AVT*^^) as a function of Oi for n = 3, = — 1 and T = 
0.155,0.165,0.175,0.185. 



5 A first look at transport properties 

One of the most interesting apphcations of holographic techniques to strongly coupled 
field theories is the computation of transport properties. In the presence of an external, 
time-dependent electric field of the form ~ e~*'^*, the conductivity reads (see e.g [HI [7]) 

a[ijj) = . (36) 
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We can use the holographic dual to compute the Fourier transform of the retarded Green's 
function for our CFT^-i at zero spatial momentum. Following the usual AdS/CFT dictio- 
nary, this boils down to computing the fluctuation of the now time dependent compo- 
nent of the bulk gauge field in the condensed background (i.e with a non-trivial profile of 
The boundary behavior of such a fluctuation is 



4(1) 



(37) 



From here, we can determine the Green's function, which results in the conductivity 



i Ax 



(3^ 



In summary, the main task is to solve the equations of motion for the A^ fluctuation in 
the condensed background. 



+ y 4 + ( ^2 f 



UJ 



2^' 



Ax = 0, 



(39) 



where the prime denotes derivative with respect to r. 

In addition, we have to impose causal boundary conditions, i.e. incoming boundary 
conditions at the horizon. This requires that, close to the horizon. 



Ax ~ 



1 



(40) 



The conductivity for n = 2 was worked out in great detail in [H] so we will focus on 
n = 3. In figure 12, we plot the real part of the conductivity as a function of u for both 



Oi and O2 condensates at various temperatures. 




<0i> 



Figure 12: Re [a] as a function of uj and T. 
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Figure 12 clearly shows the existence of a gap. The most interesting difference with 
respect to n = 2 is the appearance of resonances at certain non zero frequencies. These 
additional poles suggest the existence of bound states caused by strongly interacting low 
energy excitations (quasiparticles) of this superconductor. This behavior is qualitatively 
similar to the one observed in |27j for holographic superconductors in various dimensions 
as the bulk mass of the scalar approaches the BF bound. It would be interesting to analyze 
the transport properties of this model in further detail. 



6 More general models 

In this section we explore models with more general forms of the function JF. For this 
purpose, let us focus on 

T=^^ + C3^=^ + Ci^\ (41) 

In order to avoid keeping track of the absolute value in (|3]), we will consider C4 > and 
values of C3 such that JF > for all 



Figure 13 a shows the condensate as a function of temperature for 03,04 > 0. The free 



energy is plotted in Figure 13 b. We conclude there is a first order phase transition with 
a positive condensate. The non-vanishing C3 breaks the \E' — \1/ symmetry of previous 
examples. As we discussed around (|3]), it is sometimes possible to constraint > by 
deriving the model from one with a complex field. In this case, only the solid pieces of the 



curves in Figure 13 Complex models with the desired properties indeed exist for (41); we 



do not present them here for simplicity. 




0.29 0.30 
T 



(a) 



Figure 13: (a) Oi condensate as a function of temperature at fixed /i = — 1 for C3 = — 1 
and C4 = 1/4. (b) Normalized Aiy^^^^ at fixed ji = —1 as a function of the condensate for 
T = 0.270, 0.280, 0.305, 0.314 and 0.330. 



Let us now consider negative C3. Figure 14 a is similar to the corresponding one for 



positive C3, but shifted towards negative Oi. Figure [T4|b shows the free energy. If \1/ can 
take any real value, we conclude there is a first order phase transition with a negative 
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condensate. As before, if \1/ > only the solid branches survive. Interestingly, in this case 
the transition is second orderF°] 




0.28 0.29 0.30 0.31 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

T <0i> 



(a) (6) 

Figure 14: (a) Oi condensate as a function of temperature at fixed yU = — 1 for C3 = — 1 
and C4 = 1/4. (b) Normalized AVT*-^^ at fixed /x = —1 as a function of the condensate for 
T = 0.270, 0.280, 0.305, 0.314 and 0.330. 



In the case of second order phase transitions, it is of great interest to determine the 
critical exponent (3 that controls the vanishing of the condensate as the temperature ap- 
proaches Tc 

0,(T)~A(T,-T)^. (42) 

The n = 2 model was analyzed in where it was found that (3 = 1/2 for both Oi and 
(92j3 This value of the critical exponent is the same one that appears in the mean field 
approximation of Landau-Ginzburg. In the present case, the condensate vanishes as 

Oi(T)~0;(T,)(T,-T), (43) 

since 0[(Tc), the derivative of the condensate with respect to the temperature evaluated 
at Tc, is finite. Namely, this corresponds to f3 = 1. This behavior just follows from 
vertically shifting the condensate curve and restricting it to the positive branch. It would 
be interesting to investigate whether it is possible to find more tunable ways of obtaining 
non mean-field critical exponents. 



7 Conclusions 



We have introduced a general class of strongly coupled CFTd_i's with dual holographic 
descriptions which exhibit spontaneous symmetry breaking of a global U{1) symmetry at 

^"Notice that Oi = is a critical point of the free energy for both positive and negative C3, with no 
apparent non-analyticity when keeping the positive branch. 

^-'^Other critical exponents of this model were studied in PSI [^ . 
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low temperatures. This class contains and generalizes the model in [9]. Our dual models 
are defined in terms of a Stiickelberg-like lagrangian coupled to gravity which is determined 
by a function T . We considered some basic forms of JF, given by order n monomials and 
combinations of them. 

Interestingly, for n > 3 the phase transition is of first order, while for n = 2 it is of 
second order. We have explicitly seen this by studying the free energy, whose derivative 
exhibits the characteristic discontinuity for > 3. Interestingly, met ast ability, one of the 
characteristic features of a first order phase transition, is also clearly observed. 

We have initiated the study of the transport properties of these systems. The most 
striking feature is the appearance of extra resonances. Perhaps this feature is not so 
surprising, since for n > 3 we can understand the "vertex" '^'^[d'p — Af' along the lines 
of [30] as leading to inelastic scattering. Therefore, it is plausible that for n > 3 the 
condensate, represented by its dual field does actually have some internal structure 
which in particular allows its breaking. This would show up in the conductivity as extra 
poles. It would be interesting to understand the dependence of the gap on n for monomial 
JF or a more detailed calculation of the conductivity in order to determine whether it 
reaches some limiting curve. 

We have taken a step towards the study of general JF, by considering the model in (41 ). 
Interestingly, when the condensate is constrained to be positive because the model comes 
from one with a complex scalar, it has a second order phase transition with a critical 
exponent (3 = 1. This behavior is different from the 1/2 exponent characteristic of the 
mean field approximation that arises for the pure n = 2 model. It would certainly be 
interesting to study whether, either using models with generalized kinetic terms as the 
ones introduced in this paper or by other means, it is possible to obtain other critical 
exponents. Finally, one of the most exciting remaining questions is how the details of the 
function JF are mapped to the dual CFT. 
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